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ABSTRACT 



We analyse high-resolution near-UV and optical spectra of the afterglow of GRB 080310, obtained with the Very Large 
Telescope Ultraviolet and Visual Echelle Spectrograph (VLT/UVES), to investigate the circumburst environment and 
_ . . the interstellar medium of the gamma-ray burst (GRB) host galaxy. The VLT rapid-response mode (RRM) enabled the 

V,^ ■ observations to start only 13 minutes after the Swift trigger and a series of four exposures to be collected before dawn. 

A low neutral-hydrogen column-density (logA(Hl) = 18.7) is measured at the host-galaxy redshift of z = 2.42743. At 
this redshift, we also detect a large number of resonance ground-state absorption lines (e.g., Cll, Mgll, Alll, Sill, Cm, 
Civ, Siiv), as well as time- varying absorption from the fine-structure levels of Fen. Time- varying absorption from a 
Q ■ highly excited Fe ill energy level (^83), giving rise to the so-called UV 34 line triplet, is also detected, for the first time in 

' a GRB afterglow. The Cm ground-state and all observed Fell energy levels are found to depopulate with time, whilst 

tj^ , the Felll ^Ss level is increasingly populated. This absorption-line variability is clear evidence of ionization by the GRB, 

d ■ which is for the first time conclusively observed in a GRB afterglow spectrum. We derive ionic column densities at each 

epoch of observations by fitting absorption lines with a four-component Voigt-profile model. We perform CLOUDY 
photo-ionization modelling of the expected pre-burst ionic column densities, to estimate that, before the onset of the 
. burst, [C/H] = -1.3 ± 0.2, [0/H] < -0.8, [Si/H] = -1.2 ± 0.2, [Cr/H] = +0.7 ± 0.2, and [Fe/H] = +0.2 + 0.2 for the 

1^ . integrated line profile, indicating strong overabundances of iron and chromium. For one of the components, we observe 

vN ■ even more extreme ratios of [Si/Fe] < —1.47 and [C/Fe] < —1.74. These peculiar chemical abundances cannot easily be 

explained by current models of supernova yields. They are indicative of a low dust content, whilst dust destruction may 
also contribute to the marked Fe and Cr overabundances. The overall high iron enhancement along the line-of-sight 
suggests that there has been negligible recent star formation in the host galaxy. Thus, the occurrence of a GRB indicates 



■ that there has been episodes of massive star formation in the GRB region. 



Key words. Gamma-ray burst: individual: GRB 080310 - Galaxies: ISM - Galaxies: abundances - Galaxies: quasars: 
absorption lines 

1. Introduction In particular, optical and ncar-ultraviolct (near-UV) spec- 

troscopy of long-duration {tohs > 2 s) GRBs reveals the 



X: 

^ ■ Gamma-ray burst (GRB) afterglows, which are produced chemical composition of the"°gas along the line-of-sight 
- ■ ■ by the shock between the burst jet and either the sur- through ab sorption lines imp rinted in the afterglow spec- 
rounding intersteUar medium (ISM) or stellar wind, ra- ^^um (e.g.. IProchaska et aWmOlh . However, the vast ma- 
diate their extremely powerful and featureless continuum jprity of spectra have to date been taken at low resolution 
through their host galaxies (see |Piran| 1200J; i Meszaroj (e.g.. iFvnlio 6^111 [20091: IChristensen et al.ll2011h . hamper- 
I2OO6I for reviews). This offers a unique opportunity to j-^g nietallicity studies 
study the properties of such distant and faint galaxies. 

At the high-resolution end {R ^ 45,000), the Very 

^ „ ; ~ ~7"^ . , . . , Large Telescope (VL T) Ultraviolet and Visual Echelle 

Based on Target-OLOpportumty observations carried out Spectrograph (UVES. iDekkrel^ [2000h an d the Keck 

m service mode using the Very Large telescope Rapid-Response , „ i . 1 ', ,^ 177 — ; — ; — rn 

Mode under programme ID 080.D-0526, P.I. Vreeswijk, with the High-Resolution Echelle Spectrometer (HIRES, I Vogt et al. | 

Ultraviolet and Visual Echelle Spectrograph (UVES) installed UMM allow detailed studies of the metal abundances, 

at the Nasmyth-B focus of the VLT, Unit 2, Kueyen, operated physical states, and gas kinematics within the host 

by the European Southern Observatory (ESO) on Cerro Paranal galaxies and poss i bly lower-redshif t inte rvening systems 

in Chile. (|Fiore et all 120051: IProchaska et al.l 120071: IVreeswiik erall 
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2007t iFox et all l2008t iLedoux et all 120091: IVergani et all 
2009t) . albeit for a small sample of bright afterglows. Among 
the six Swift-cia. GRB afterglows with metallicity deter- 
minations from VLT/UVES spectroscopy, four absorbers 
have low metallicities, [X/H] < —1, where X is taken to 
be either Zn or S. These generally have low dust deple- 
tion factors, [X/Fe], indicating a low dust content, while in 
one case a hig her [S/Fe] ratio is sug gestive of an a-element 
enhancement ([Ledoux et alj 12009^ . On the other hand, 



Table 1. Solar abundances used in this paper 



iProchaska et all (|2007^ found that [a/Fe] and [Zn/Fe] tend 
to be higher in GRB absorbers than along quasi-stellar ob- 
ject (QSO) lines-of-sight, based on high- and low-resolution 
spectroscopy of a sample of 16 absorbers, indicating signif- 
icant contributions from massive stars and/or differential 
dust depletion. However, the current GRB absorber sam- 
ples are still too limited in size and possibly too biased to 
determine what the typical abundances of GRB host galax- 
ies are. 

A larger collection of metal column densities in 29 GRB 
absorbers , drawn from both high- and low- resolution spec- 
troscopy (jSchadv et al.ll201]] ). reveals that the ranges of rel- 
ative chemical abundances measured in GRB absorbers are 
generally similar to those determined along QSO lines-of- 
sight. In this paper, we present peculiar metal abundances 
in the host galaxy of GRB 080310 and discuss their possible 
origin. 

Thanks to time-resolved high-resolution spectroscopy, 
the variability of fine-structure lines has been detected in 
GRB absorbers and explained by photo-excitation induced 
by the incident GRB afterglow flux (i.e., UV pumping, 
IVreeswiik eralll2007l: IProchaska et al.ll2006[) . On the other 
hand, the variability of resonance metal lines (i.e., those 
associated with the ground-s tate level) is expecte d in the 
case of photo-ionization (e.g.. lPerna fc Lo"eblll998D but has 
never been conclusively detected until now (the possible 
variability of Lya towards GRB 090426 was reported by 
[Th one "etall[201ll) . Here we present the detection of signif- 
icant variability in both fine-structure and resonance metal 
absorption lines in the afterglow spectrum of GRB 080310. 
In addition, a triplet of lines arising from a highly excited 
level of Fe iii is observed in this absorber, which is a primer 
in a GRB afterglow spectrum. As we discuss in this paper, 
these observations clearly show that photo-ionization by 
the GRB has taken place. Self-consistent photo-excitation 
and ionization modelling to determine the GRB-absorbing 
cloud distance will b e presented in a companion paper 
(jVreeswiik et al.ll2012L hereafter referred to as paper II). 

This paper is organized as follows. In Sect. 2, we re- 
port on the observations and the data reduction. In Sect. 3, 
we derive ionic column densities from Voigt-profile fitting 
of identified absorption lines and in Sect. 4 we present 
CLOUDY photo-ionization modelling. We discuss our re- 
sults in Sect. 5 and conclude in Sect. 6. Throughout the 
paper, we adopt ions cm~^ as the linear unit of column den- 
sity N. The relative abundance of two chemical elements, X 
and y, is defined as [X/^] = log — log %[y)q ■ ^® 
timate its uncertainty by adding in quadrature the errors in 
the observed column densities and in the solar abundances 
involved. For the reference solar abundances appearing in 
the secon d term of this formula , we follow the recommen- 
dations of iLodders et al.l ()2009f ). adopting e ither their me- 
teoritic estimates, the photospheric values of lAsplund et al.l 



Element 
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S" 


Element 


AiE\) ± a A " 
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Cr 


5.64 ± 0.04 
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7.47 ± 0.04 
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8.69 ± 0.05 
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Ni 
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7.51 ±0.01 
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4.63 ± 0.04 
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Not es. " Abundanc e s A (E l) = log N(E1)/N(H) ± 12 are taken 
from lAsplund et al.l (|2009D following the recommendations of 
ILodders et al.l (|2009r i. Source of the estimate: solar photo- 
sphere (s), meteorites (m), or the average between the two (a). 



(|2009D . or the average of these both. The solar abundances 
used in this paper are summarized in Table [1] 

2. Observations and data reduction 

GRB 080310 triggered the Swift Burst Alert Telescope 
(BAT) on March 10 2008 at UT 08:37:58 (jCummings et al.l 
'20081). One and a half minutes later, the Swift X-Ray 
Telescope (XRT) provided an X-ray afterglow position of 
RA = 14h40'"13.5" and Dec = -00°10'32.1" (J 2000), with 
an uncertainty of 5". This detection led us to trigger the 
RRM of VLT/UVES, as set up by our team, to slew the 
Kueyen telescope to the target automatically. Once the ob- 
ject was centred on the slit, the UVES observations started 
(UT 08:51:01 ), i.e., 13 minutes and three seconds after the 
BAT trigger (jVreeswiik et al.ll2008D . 

Four exposures were taken, of increasing durations to 
compensate for the fading afterglow brightness, alternat- 
ing between dichroic #1 and dichroic #2 with standard 
UVES settings 346 + 580 and 437 + 860, respectively. In 
this way, the near-UV/optical range was basically entirely 
covered via two consecutive exposures. The instrument set- 
tings used have three overlapping regions of widths ^ 150 
A, ~ 230 A, and ^ 240 A, which were observed at all 
four epochs, and two '-^ 100 A gaps in the red owing 
to the physical separation between the two red CCDs. 
The CCDs were set to bin pixels 2x2 and the spec- 
trograph slits to be 1" wide, providing a mean resolving 
power R = A/AA ~ 46,400 (full width at half maximum 
FWHM~ 6.5 km s-i). 

The data were reduced using the ESO UVES pipeline 
v2.9.7 based on MIDAS (jBallester et al.l |2000D and the 
wavelength scale was then transformed to the heliocentric 
rest-frame. As a consistency check on the quality of the 
reduced spectra, an independent data reduction was also 
performed using the UVES pipeline v3.4.5 based on the 
Common Pipeline Library (CPL)13, providing very similar 
quality products. The observing log and details on the ob- 
servations are presented in Table |5J 



3. Absorption-line analysis 

3.1. Detected lines 

We first combined the spectra obtained at the four differ- 
ent epochs of observations, to achieve the highest signal- 



^ ESO MIDAS and CPL- based UVE S pipelines are available 
at ,http: / /www.eso.org/ sci/data-processing/ software/pipelinesT] 
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Table 2. Log of VLT/UVES observations 



Ep. 


,a 

••start 


At'' 


^cxp 


Die. 


Setting 




Coverage 




FWHM'^ 


Mean 


S/N'' 






UT 


(min) 


(min) 


# 


Ac (nm) 




Aobs 


(nm) 




blue/red (") 


airmass 


blue/red 


blue/red 


I 


08:51:05 


14.55 


3 


1 


346+580 


303- 


-388; 


476- 


684 


1.1/1.0 


1.12 


1.9-2.0/4.0-4.5 


48,500/45,600 


II 


08:56:46 


21.23 


5 


2 


4374-860 


373- 


-499; 


660- 


1060 


1.2/0.9 


1.12 


3.4-4.7/2.7-8.4 


47,900/45,300 


III 


09:04:20 


31.25 


10 


1 


3464-580 


303^ 


-388; 


476- 


684 


1.5/1.3 


1.14 


3.1-3.7/6.8-7.6 


48,500/45,600 


IV 


09:16:36 


50.06 


23 


2 


437-1-860 


373- 


-499; 


660- 


1060 


1.6/1.2 


1.17 


6.0-8.7/4.8-14.3 


47,900/45,300 



Notes. " March 10 2008. Mid -exposure time after the burst event in the observer's frame. Full width at half maximum of 
the spatial profile in the two-dimensional spectrum. Signal-to-noise per pixel; min-max ranges over the spectral regions where 
absorption lines are observed. 



Table 3. Transitions simultaneously modelled with a Voigt 
profile. 



Epochs 


Ion 


Level 


Transitions A (A) 




Sin 


2p 

-T 1/2 


1526 




Fell 


D9/2 


1608 


I, III 


Felll 


'D4 


1122 




Felll^'^*" 


^S3 


UV34: 1895, 1914, 1926 




Cll 


'Pl/2 


1334 




Sin 


2p 
-t^l/2 


1190, 1193, 1260, 1304 




Sin* 


^P3/2 


1264, 1265 




Fen 


D9/2 


2344, 2374, 2382*, 2586 


II, IV 


Fen* 


D7/2 


2333, 2365, 2389, 2396, 2612 




Fen** 




2349, 2399, 2405 




Fen*** 




2338, 2407, 2411 








2411, 2414 




Fcni 


'D4 


1122 



Notes. " Seventeenth energy level above the ground-state. ' 
Only the redder component is fitted since the absorption profile 
is saturated in bluer regions. 



to-noise ratio (S/N) for the identification of absorption fea- 
tures. Starting with the detections of the strongest dou- 
blet lines and confirming these with associated resonance 
lines detected at the same redshift, we identified absorp- 
tion lines belonging to eight distinct systems. We associated 
the highest redshift system, at z = 2.42743, with the GRB 
host galaxy. The presence of Fe 11 fine-structure lines in this 
system confirms this association , since their origin is m ost 
likely linked to the GRB (e.g., Ivreeswiik et al.ll2007t) . In 
this paper, we focus our analysis on the absorption within 
the host galaxy, and exclude the absorption lines associated 
with the seven intervening systems. The latter are listed in 
Table lA.ll (see Appendix) with equivalent width (W) esti- 
mates. 

At the GRB redshift, we confirm the presence of absorp- 
tion features such as Lya, and highly ioni zed specie s , such 
as C iv, Siiv, N V and O vi, as reported in lFox et all (|2008[ ) 
and lLedoux et al.l()2009D . We also detect low- ionization ab- 
sorption lines, such as Cii A 1334, Alii A 1670, Aim AA 
1854, 1862, Mgii AA 2796, 2803, and a blended Oi A 1302, 
in addition to Si 11, Fe 11 ground-state and fine-structure, 
and Fe ill ground-state transitions. Moreover, three lines 



displaying similar profiles are clearly detected, but cannot 
be associated with any transition reported for GRB lines- 
of-sight to date, nor any of the intervening systems. We 
identify these transitions as the Feiii UV34 triplet AAA 
1895, 1914, 1926. We discuss these lines and their excited- 
level origin in detail in Sect.O Remarkably, the line profiles 
at different epochs reveal the clear variability in both the 
Fell and Feiii transitions. Fig. [T] shows the (non-) vari- 
ability in some representative transitions, in the left panel, 
and the Fe 11 and Fe ill apparent column densities (derived 
by combining together the transitions of each level), in the 
right panel. A decrease in the Fe 11 apparent column den- 
sity is evident in all the levels, while the excited Feiii ""Ss 
population increases with time. We discuss the evidence for 
the variability of Fe 11 and Fe ill and other species further 
in Sect. |331 

3.2. Voigt- profile fitting 

Table [3] lists all metal lines, including both ground-state 
and excited-level transitions, that are clearly detected yet 
not strongly saturated and therefore suitable for the Voigt- 
profile fit. A first look at the absorption lines identified in 
the spectrum reveals that their profiles are not identical 
for all the ions. In particular, the strongest component of 
all the Fe 11 lines appears broader and blue-shifted by ^ 20 
km s~^ with respect to the strongest component of the line 
profile of the other low-ionization species, such as Si 11, C 11, 
and Aim (see, e.g.. Fig. [T|). 

Before investigating the physical origin of this shift, we 
need to ensure it is not an artefact caused by the wave- 
length calibration. In this check, a first point of reference is 
the Fe iii A 1122 transition, which is observed with a similar 
profile at the same wavelength, despite the noise, in all four 
spectra in the overlapping region between the two different 
settings. Second, while most of the Fell transitions are ob- 
served in the red arm spectra obtained with dichroic #2, 
the bluer Fell A 1608 shows the same Fell profile in the 
spectra taken with dichroic #1. Moreover, the Fell A 1608 
and the Si 11 A 1526 lines are observed simultaneously within 
the same spectral region, indicating that the wavelength dif- 
ference between the strongest iron and silicon components 
is real. In addition, the Fell AA 1144, 1096 transitions that 
are observed in the blue arm spectra obtained with the 
dichroic #1 and #2 show the same kinematics as the red 
arm Fell (stronger component at —20 km s~^). In particu- 
lar, the Fell A 1144 shows clear variability (see Fig.[T]), con- 
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Fig. 1. Left: a selection of line profiles observed at different epochs (epochs I and II in blue; epochs III and IV in black). 
Telluric features are highlighted in yellow. The dashed lines mark the position of the strongest component of all the Fe II 
lines and the strongest component of Sill, Cii and Aim lines, respectively. Right: the apparent column-density profiles 
(smoothed with a boxcar of five pixels) of the Fell and Fciii levels, obtained by combining together different transitions 
that probe a particular level, show a clear evolution with time. In particular, while all Fell levels depopulate, Feiii shows 
the opposite trend. 



firming the line identification. Finally, a synthetic telluric 
spectrum (see next Section) closely matches the observa- 
tions across the full spectral range, proving the reliability 
of our wavelength calibration. 



The absorption-line profiles of these low-ionization 
species diffe r significantly fro m the high-ionization ones 
analysed bv iFox et"all (|2008[ ). with Civ, Ovi, and Siiv 
having much broader line profiles with blue-shifted wings 
(up to 220 km s~^ for Civ). Thus, we can exclude a di- 
rect association between the high-ionization wings and the 
low-ionization components. A physical separation between 
the low- and high-ionization species, and therefore different 
line profiles, is not unexpected, as they may trace different 
regions of the ISM owing to their different ionization po- 
tentials. 



3.2.1. Cii, Aim, Sin, Feii, and Feiii 

To derive the column densities of the metals and their evo- 
lution with time, we modelled the observed lines with the 
Voigt-profile fitting software VPFITQ. We first normalized 
the spectrum locally using carefully selected featureless re- 
gions, around each line, by fitting a zero- or first-order poly- 
nomial to these regions. We considered the best-fit contin- 
uum level plus and minus 0.5 times the noise root mean 
square in the ad jacent continuum (u pper and lower con- 
tinua), following iLedoux et al.l (l2009l ). This allowed us to 
take into account the normalization error in the estimate 
of the error in the column densities. The normalization er- 
ror tends to dominate over the formal Voigt-profile fit error 
for weak transitions. 

To avoid false identifications caused by sky absorption 
lines, we created synthetic telluric absorption-line spec- 
tra (see, e.g.. Fig. [T]). These telluric spectra were obtained 



http : //www. ast . cam. ac .uk/'^rf c/vpf it .html 
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with an IDL routine acting as a wrapper to the Reference 
Forward Model (RFM) version 4.28, a hne-by-line radiative 
transfer cod^ The HITRAN'2008 (|Rothman et al.ll2009[ ) 
database provided the molecular absorption-line parame- 
ters. The atmospheric profiles describing the pressure, tem- 
perature, and humidity at 20 different atmospheric layers 
were obtained from the Air Resources Labora tory READY 
Archived Meteorology websitcQ. Following ISmette et al.l 
()201Cll) . we calculated synthetic telluric absorption-line 
spectra for different amounts of precipitable water vapour. 
Fitting the synthetic telluric spectra to the observed spec- 
tra between 7100 A and 7450 A provides a best-fit pre- 
cipitable water-vapour column of 3.84 ± 0.05 mm. In the 
Voigt-profile fits, we systematically discarded the spectral 
regions where the telluric absorption features drop below 
99% of the continuum level. 

We treated the spectra taken at each epoch separately, 
to allow for the detection of absorption-line variability. 
Most species are observed only at two epochs, since the 
wavelength coverage differs between epochs. Exceptions are 
the Feiii A 1122 resonance line and its fine structure Fciii* 
A 11243, which are covered at all four epochs. To simultane- 
ously fit all the line profiles listed in Table [3] with the same 
components, we first focused on the highest S/N (epoch IV) 
spectrum. We then included the epoch-Ill spectrum to also 
cover the Feiii UV34 triplet. 

Noticeably, the lines associated with different ions show 
different profiles, making their decomposition non-trivial. 
As a first approach, only the stronger Cil, Sill, Sin*, and 
Fell ground-state lines were considered for the identifica- 
tion of the four velocity components "a" , "b" , "c" , and "d" , 
listed in Tabled that are necessary to describe the line pro- 
file adequately. Below, we briefly discuss the justiflcation for 
this decomposition. 

First, the C ii profile is mainly described by a strong 
narrow (b = 7 km s~^) component, which we refer to as 
component "c", at redshift z = 2.42743 ± 0.00002 0. We 
adopt this redshift as our zero- velocity reference point. At 
this same position, we detect the strongest component of 
the Sin lines, as well as all the other identified transitions. 
However, despite the strength of component "c" in terms 
of silicon and carbon, iron does not show its stronger ab- 
sorption at this wavelength position. Second, the Fe n main 
contribution, which we refer to as component "b" , is blue- 
shifted with respect to the main Si and C components, but 
is also clearly observed in C and Si. Third, the Fe n profiles 
display an even bluer, quite- weak-but-evident, absorption- 
component "a" , that is almost undetected or very weak for 
the other ions. Finally, a redder component, "d" , widens 
the profile of the Si lines, although this is less obvious for 
the Fen lines. To investigate the contribution of compo- 
nent "d" to Fen we include in the fit only the red wing 
of Fen A2382, which is otherwise saturated (Fen A2600 is 
also saturated). Given the resolution of the instrument of 

^ http://www.atm.ox.ac.uk/RFM/ 

http : //ready . arl .noaa.gov/READYinf o .php 

^ We use the non-detection of Felll* A 1124 to put an upper 
limit on the Felll* population. 

^ The dominating systematic error, arising from the heliocen- 
tric velocity and wavelength calibration, is estimated by i) com- 
paring our spectra with a synthetic telluric one and ii) mod- 
elling simultaneously several transitions from different regions 
of the spectrum, observed with different dichroics and at differ- 
ent epochs. 



^ 7 km s~^, the individual components derived with the 
Voigt-profile fit are blended, and therefore the decomposi- 
tion is not unique. However, the remarkable difference be- 
tween the Fe n transitions and the Si n and C n lines enables 
component "b" and component "c" to be distinguished and 
therefore more reliably constrained. 

We then fixed the positions (rcdshifts) and widths 
(velocity-broadening parameter) of components "a","b", 
"c" , and "d" , reported in Table|31 and included all the tran- 
sitions listed in Table [3] at their highest S/N (all transitions 
detected at epoch IV and the Fe ni transitions at epoch HI) , 
in a simultaneous four-component fit of all lines. A pure 
turbulent broadening (i.e., neglecting thermal broadening) 
provided the best fit to the line profiles. Figs. [2] and [3] show 
the model fit, along with its individual components, con- 
volved with the instrumental resolution; the reduced chi- 
square is xt = 1-4- As mentioned above, the spectral re- 
gions affected by telluric features were excluded from the 
fit. 

To combine the errors derived from the Voigt-profile fit 
with the errors introduced by the normalization, we mea- 
sured the column densities for the three continuum levels 
(best fit, lower, and upper) and we add in quadrature the 
difference between the column densities of the various con- 
tinua with the formal VPFIT errors. The four-component 
decomposition is significant for all the transitions in the 
upper continuum level, i.e. all components are required by 
the software to describe the line profiles. In the case of 
a non-detected component in the best-fit continuum, and 
for weak components with upper errors of more than twice 
the inferred column density, we used an upper limit corre- 
sponding to the largest value between the column density 
estimate from the upper continuum level and the 3cr de- 
tection limit. For components that are well-constrained in 
the best-fit continuum, but not detected in the lower con- 
tinuum, we adopted the expectation value of half the 3tT 
detection limit to calculate the lower error. The detection 
limits for each component were calculated by measuring 
the uncertainty in the W from the error spectrurrQ at the 
position of the stronger transition of each ion, and con- 
verting it to column densities assuming the optically thin 
limit. An aperture of size twice the FWHM was chosen to 
enclose ~ 98% of the surface, for each component found in 
the Voigt-profile fit. We used FWHM, = 2^2 ln2(Ti, where 
erf = bf/2 is the Gaussian variance and bi is the Doppler- 
broadening parameter of component i. Since the proxim- 
ity of components "c" and "d" may cause some degener- 
acy in the measurement of the individual column densities, 
we also considered the results of the combined components 
"c-|-d" in our abundance analysis. The 3(t detection limit 
to components "c-|-d" and the total line profile are mea- 
sured assuming the overall apertures necessary to cover the 
blended components, i.e., two FWHMs correspond to 27 
km s"""^ and 110 km s~^ for "c-|-d" and the total profile, 
respectively, corresponding to 6 = 16 km s~^ and 66 km 

We applied the model fit derived above to the spec- 
tra collected at all the other epochs, allowing only the col- 
umn density to change, in order to investigate any spectral 
variability. As a consistency check, an independent fit of 



^ AWrest = pixel size [A] x \/Ei normalized error^ /(I -I- z), 
where i spans the chosen aperture. 
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Fig. 2. Voigt-profile fit of absorption lines observed at epocli IV. Tlie black solid lines show the observed spectrum, 
with the associated error indicated by the dashed line at the bottom. The red lines show the combined model fit. Other 
coloured lines show the individual components. The four components, "a", "b", "c", and "d", are labelled accordingly. 
J'he transitions are listed in order of decreasing oscillator strength, for each population. Telluric features are highlighted 
in yellow. 
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Table 4. Column densities 
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Notes. " We adopt z = 2.42743 as our zero-velocity reference point. Purely turbulent broadening. Logarithm of the 
ion column density and la error estimate or 3cr upper limit. The detection criteria for each component are: i) detection in 
at least both best-fit and upper continua; ii) the column density is larger than the 3a detection limit; iii) its total error 

Square brackets indicate the transitions used to measure the 3a upper lim- 
its. The number of *'s indicates the excited energy level above the ground-state, either fine-structures of the ground-term or 
higher excited levels. ^ Conservative limits to account for any possible underestimate of Sin and Cll in component "c" in case of 
a low 6- value. ® Upper limit derived from deblending, see Sect. 13.231 



errtot \/ ^'^'"formal ^^^normalization ^'-'S2. 
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Fig. 3. Voigt-profilc fit of absorption lines observed at epoch. III. Tlie fit decomposition was derived from tire highest 
S/N spectrum (obtained at epoch IV). 



epoch-II data confirms the positions and the widths of the 
components. 

Tabled lists the column densities derived from this four- 
component fit, for the four epochs. The summed compo- 
nents "c-|-d" and the sum of all the components are listed 
in the last two columns, for each level. We conservatively 
considered lower limits to the C ii content in component "c" 
because of the risk of the C ii A1334 line being saturated 
with a 6-value possibly lower than 7 km for this com- 
ponent. However, the non-detection of Cii* A1335 - com- 
monly detected in QSO damped Lyg systems (QSO-DLAs, 
logA^(Hi) > 20.3. IWolfe et al.lll986[ ) - also indicates a low 
carbon content. In the same way, we conservatively consid- 
ered lower limits to the Si ii column densities of component 
"c" at epochs II and IV, when the stronger Sin A1260 line 
was included. However, the Si ii contents estimated from 
epoch-I and cpoch-III data do not differ much from those 
of epochs II and IV, suggesting that any Si ii underestimate 
must be small. 

When summing upper limits, we used the non- 
significant estimates from the individual components when 
available!!, and otherwise expectation values of half the up- 
per limits. This takes into account that individual compo- 
nents could be non-significant when considered singularly 
but that their sum may be 3(t significant. If none of the in- 
dividual components to be combined have a non-significant 
column-density estimate, then the 3tT detection limit to ei- 
ther "c-|-d" or the total line profile is used. As noticed be- 
fore, some transitions are covered at only two epochs. 

In addition, we summed the column densities of all the 
levels (i.e., ground-state and excited levels) belonging to 
the same ion (Ol, Sin, Fen, and Fcni) and list them in 

* We only consider those non-significant individual- 
component column- density estimates that have an error 
smaller than aiogjv < 1, to avoid including measurements that 
are totally unconstrained. 



Table 5. Total ionic column densities including excited 
states. 



Ep. 


Ion 


comp. a 


logA^i o-r„gjv 
comp. b comps. c+d 


total 


I 
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< 14.46 


14.081°;^ 
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II 
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Sin 
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III 
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< 14.04 
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IV 


Cii 
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Sin 
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< 13.23 
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13.59trS 
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< 14.17'' 

12.66tO;}2 
n 97+0 07 


> 14.12 
< 13.84 

> 13.41 
13.50l°:»? 


> 14.17 
< 14.39 

> 13.51 

14 20+" °^ 



Notes. " Logarithm of the ion column density and correspond- 
ing la error estimate. ^''' Upper limits from deblending. The 
highlighted column stresses the results for component "b" , which 
shows the most peculiar [Sill/Fell]. 



Table [51 We again present the results for the individual 
and combined components. 



3.2.2. Crn 

In addition to the low-ionization lines identified in the pre- 
vious subsections, we identified Crii AA2056,2062 absorp- 
tion features at —20 km s'^, corresponding to component 
"b" in the line profile modelled for the other metals. We 
fitted the Voigt profile of component "b" (fixed z and b- 
value) to these two transitions, together with the weaker 
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Fig. 4. The Cm AAA2056,2062,2066 transitions at epochs 
II (blue) and IV (black). The component "b" Voigt-profile 
fit at epoch II is shown by the solid red curve. The er- 
ror spectra are displayed at the bottom of each panel with 
dashed lines. 



Cm A2066, deriving the column densities listed in Table SI 
The other three components were not detected in the spec- 
tra and therefore we calculated the upper limit to the col- 
umn densities as described above for the other transitions. 
We treated the Cr transitions separately from the other ele- 
ments in order not to bias our results, as they are weak and 
only detected in component "b" . Remarkably, the Cr ii col- 
umn density for this component varied with time, at a 2.2 a 
level, as calculated from the Cr ii column density at epoch 
II and the expectation value on the upper limit at epoch 
IV. However, this is a conservative approach, as the errors 
in the column densities include the errors on the normaliza- 
tion. Using only the formal errors from the Voigt-profile fit, 
we derived a 4.6 a decrease in the Cr ii column density from 
epoch II to epoch IV (logiVepii(Crli) = 13.07 ± 0.06 and 
logA'^Bpiv(Cm) = 12.37 ± 0.14). This is the first time that 
Cr II absorption has been observed to vary in a GRB host 
galaxy. The spectral variability of Cm AAA2056,2062,2066 
is shown in Fig. 21 

3.2.3. Hi and Oi 

We derived the neutral-hydrogen column-density log A^(II i) 
= 18.7 ±0.2 and logA^(Hl) = 18.7 ±0.1 at epochs II and 
IV, respectively, from a Voigt-profile fit of the Lya ab- 
sorp tion line (Fig. [3, using the MIDAS FITLYMAN pack- 
age (jFontana fc Ballesterjll995[ ). These values were derived 
from the Voigt-profile fit of the red wing of Lya, fixing 
the central redshift and the 6-value to the one found for 
the metal absorption system (mean of components "b" and 
"c" ) . The extra absorption in the blue wing of the line could 
be related to either the GR B host-galaxy h alo or galaxy 
winds observed in Civ (see iFox et al.l[2b08^ . As a sanity 
check, we also independently fitted the whole Lya profile 
with the redshift and 6-value as free parameters. In this 
case, a larger b was constrained by the extended core of 
the line, while TV (Hi) is constrained by the wings of the 
line profile, providing logA^(Hi) = 18.8. Thus, there can- 
not be significantly more Af(Hi) in this system than what 
is determined from the one-component fit at the redshift of 
the metal lines. This low column density classifies the ab- 




-500 500 

Relative velocity (km s~i) 

Fig. 5. A neutral- hydrogen column-density fit to the 
damped Lya line at the host-galaxy redshift, shown with 
the red solid curves, provides logAf(Hi) = 18.7 ± 0.2 at 
epoch II (top panel) and log (Hi) = 18.7 ± 0.1 at epoch 
IV (bottom panel) . The redshift was fixed to the mean red- 
shift of the low-ionization metal components "b" and "c" 
{z = 2.427315, taken here as the zero-point of the velocity 
scale) , and the velocity broadening was fixed to 6 = 20 km 
s~^, corresponding to the difference in velocities of these 
two components or alternatively the quadratic sum of their 
6- values. The la errors are shown with dashed lines. To 
show that our measurement is robust, we also performed a 
fit - shown by the green curve - in which we allowed both 
the redshift and velocity broadening to vary as free param- 
eters, finding z = 2.42670, 6 = 46 km s~\ and logiV(Hi) 
= 18.8 for both epochs. The small vertical lines mark the 
central redshifts of the two independent fits. 



sorber as a Lyman-limit system (LLS, ISargent et al.l[l"989l : 
iRauch 1998: Pcroux e t al. 2 00l with 17.3 < logAf(Hi) 
< 19. Finally, we note that the absorption feature blue- 
shifted by roughly 600 km from Lya is probably due to 
Hi absorption along the line-of-sight and possibly related 
to the host galaxy. However, this component is unrelated 
to the metal absorption system that we focus on in the 
analysis. 

The Oi A 1302 absorption associated with the GRB 
host galaxy is blended with the Alii A1670 line from the 
z = 1.6711 intervening system, listed in Table lA.ll The 
presence of Al ill absorption from the intervening system 
allows the Alii to be deblended from the GRB host-galaxy 
Oi. As the detected Aim AA1854,1862 doublet is strong, 
and Fell and Sin are also detected at this redshift, a sig- 
nificant contribution from Al II is indeed expected. A total 
ionization of Al ii can be excluded, its potential being higher 
than Fell and Sin. 

At the z = 2.42743 location of O i A1302, we observed a 
very weak absorption feature that may be associated with 
component "c" modelled as described above. However, the 
detection is formally not significant, given the noise level. 
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Table 6. Nv column densities and variability 





lo^ 


; (IN VJ± (Tiog 


JV 




blue comp. 


red comp. 


total 


Epoch II 


13.87±0.05 


13.70±0.06 


14.10±0.04 


Epoch IV 


13.81±0.02 


13.68±0.03 


14.05±0.02 


"IIJV 


1.1 


0.3 


1.1 


(km s^^) 


-53 ±1 


1 ± 1 




b (km s^^) 


23 ±2 


18 ± 2 





Notes. °' Logarithm of the N v column density and correspond- 
ing 1(7 error estimate, as derived from a Voigt-profile fit of the 
AA1238,1242 transitions {xl = 1-126 averaged between the two 
epochs). ^ Significance of the epoch II - epoch IV variation. 
Zero velocity reference aX z = 2.42743. 



No significant absorption was detected for components "a" 
and "d" either. Thus, we considered Scr detection limits to 
the O I column density for components "a" , "c" , and "d" , 
listed in Table |H To investigate how much O i from the 
host-galaxy component "b" can be hidden in the Alil ab- 
sorption associated with the intervening system, we first 
fitted the Al ii and Al ill line profiles together with a Voigt 
profile, in the upper continuum level, for epoch IV, and 
we identified two components "a" and We then in- 
cluded the host-galaxy O i contribution in the fit, resulting 
in a best-fitting xl = 1-12 for logiV(Oi, "b")= 13.50 and 
the aluminium column densities listed in Table IA.2I in the 
Appendix. The line profiles for the best fit are presented 
in Fig. IA.3I in the Appendix. For a conservative estimate 
of the Oi column-density limit, we continued to increase 
the amount of oxygen in component "b", by allowing the 
Al II decrease until the reduced xt deviated more than 3cr 
from the best fit, resulting in the upper limit logA^(Oi, 
"b")< 14.10. We repeated this procedure for epoch II. 

3.2.4. Nv 

Although IFox et ahl (|2008[ ) comprehensively discussed the 
high- ionization species along several GRB lines-of-sight, 
here we included N v because we expected it to vary, given 
the Fell and Feiii variability observed here and the ev- 
idence, in so me cases, that these tr ansitions arise in the 
GRB vicinity (jProchaska et al.ll2008D . We analysed the Nv 
associated with the GRB host galaxy independently of the 
low-ionization species, since its high-ionization state asso- 
ciates it with a different gas phase with different kinematics, 
than the low-ionization species analysed in this paper. 

The Nv AA1238,1242 doublet shows a clear two- 
component profile (see Fig. [1]), which we modelled using 
two Voigt profiles. We again used the highest S/N spec- 
trum taken at epoch IV to derive a model-fit profile to ap- 
ply to the spectra at other epochs. As a consistency check, 
we verified that a totally independent fit between epochs 
(i.e., leaving all the parameters free to vary at each epoch) 
did not change our results significantly. The Nv column 
densities at different epochs, summarized in Table IHl are 
consistent to within l.lcr, confirming the non- variation. 



Table 7. Significance level of variability in tj units 



Population 


•^var , b 


•^var, total 


Cii 


1.8 (2.8) 


0.6 


Sin 


2.1 (2.5) 


1.1 


Sin* 


0.2 (0.6) 


0.2 






9 9 


Ecu 


3.9 (3.9) 


6.3 


Eeii* 


2.5 (3.5) 


3.9 


Fen** 


3.2 (5.3) 


4.3 


Fen*** 


2.3 (5.0) 


2.8 


Fg j j ^ ^ ^ ^ 


1.1 (2.7) 


2.1 


Fein 


1.7 (1.7) 


0.6 


Eciii^'^* 


3.2 (5.4) 


2.5 



Notes. (Tvar between epoch II and epoch IV (I and III for 
Fein^^*), calculated for component "b" and for the total line 
profile, considering the conservative larger errors that include 
the normalization error (the value in parenthesis considers only 
the formal errors from the Voigt-profile fit). 

3.3. Evidence for time variability 

The Fe ii ground-state and fine-structure populations 
strongly decrease with time during the afterglow observa- 
tions, as can be seen in Fig. [1] and Table SI On the other 
hand, we observed a weak increase in the Fe ill ^83 excited 
population, that is more pronounced in component "b" , and 
in general more significant when we consider only the for- 
mal error in the Voigt-profile fit (i.e., not including the error 
in the normalization). While the Cm ground-state popu- 
lation shows hints of evolution with time, as presented in 
Sect. 13.2^ we did not find any evidence of the variability of 
Cii, Sin, or Nv (see Sect. 13. 2. 4p . In Table [71 we report the 
significance level of variability in a units, for the total line 
profile and the individual component "b" , including (and 
excluding) the error in the normalization. 



4. CLOUDY pre-burst photo-ionization modelling 

To investigate the physical conditions and metallicity of 
the host-galaxy ISM prior to the burst, w e ran a series 
of C LOUDY photo-ionization simulations (|Ferland et al.l 
Il998l ) designed to reproduce the observed ionic column den- 
sities. These time-independent (equilibrium) simulations do 
not apply once the GRB radiation has altered the ioniza- 
tion and excitation balance of the gas, but can be applied to 
the pre-burst column densities. These column densities were 
estimated from a photo-excitation/ionization modelling of 
the colum n density variab i lity o f this GRB afterglow per- 
formed in IVreeswiik et aH (j2012l paper II), to produce the 
following results log7V(Hi) = 18.7, log7V(Feii) = 14.8, 
log7V(Feiii) = 14.5, log7V(Siii) = 13.7, and log7V(Oi) 
< 14.6, which we then used as input to the CLOUDY sim- 
ulations. 

The CLOUDY models assume that the gas exists in 
a plane-parallel uniform-density slab exposed to the z = 
2.5 extragalactic backgroun d (EGB) radiation incor porated 
within CLOUDY (based on lHaardt fc MadaulfloM) . which 
has an ionizing photon density log (n.^/cm"'^) = —4.7. We 
ran a grid of simulations for different values of the ioniza- 
tion parameter U , displayed in Fig. [B] (left panel), where 
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U = n^/nH, the ratio of ionizing photon density to gas den- 
sity. We then determined (a) the best-fit value of log U by 
matching the observed Feiii/Feii ratio, whieh is a mono- 
tonically increasing function of log U (see Fig. [51 bottom 
panel); (b) the best-fit value of [Fe/H] = -1-0.2 by matching 
the Fe II and Fe ill column densities at the best-fit log U ; (c) 
the values of [C/H], [Si/H], and [Cr/H] and the upper limit 
to [0/H] by matching Cii, Sill, Cm, and Oi at the best-fit 
log U. The results are: \ogU = -3.8, [C/H] = -1.3, [0/H] 
< -0.8, [Si/H] = -1.2, and [Cr/H] = +0.7. One clear re- 
sult of this process, which we discuss in more detail below, 
is that the solar relative abundances fail to reproduce the 
column densities. We were able to estimate the error in the 
relative abundances to be ~ 0.2, given the error in the pre- 
burst [Feiii/Feii] ratio of ^ 0.1 dex (see paper H) and a 
similar uncertainty arising from the CLOUDY modelling. 

To investigate the effect of a different radiation field on 
the results of our CLOUDY model, we ran the same grid 
of simulations as described above but as suming th e Milky 
Way (MW) radiation field described in iFox et all (|2005D . 
This should account for some star formation contributing 
to the UV radiation of the host galaxy. Although GRB host 
galaxies are generally highly star-forming (with the possi- 
ble exception of GRB 080310, as discussed in Sect. [5]), their 
mass is typically ra ther low (on average M = lO^'^M©; 
ISavadio et "aD I2009D . Thus, using the integrated MW ra- 
diation field could provide a rough idea of the effect of 
star formation in the host galaxy on the ambient medium. 
The results of the alternative MW-radiation field CLOUDY 
model are shown in Fig. [5] (right panel). In particular, the 
[Fe/H] and [Cr/H] ratios are both solar and super-solar for 
the two models, varying only little with the assumed radia- 
tion field. In addition, the [Si/Fe], [C/Fe], and [Si/Cr] ratios 
are similar - peculiarly low - in the two cases, independently 
of the model assumptions. 



5. Discussion 

5.1. Low N (Hi) environment 

Most absorbers at the GRB redshift are DLAs, typically 
show i ng a higher A^(H i) than QSO-DLAs (IVreeswiik et al.l 
120041: iJakobsson et all 120061: iFvnbo et al.ll2009D . This is a 
likely consequence of GRBs occurring in star-forming re- 
gions within their hosts, while examining the spectra of 
QSOs provides information about more random lines-of - 
sight through intervening galaxies (jProchaska et al.ll2007l ). 
Only a handful of Swift GRB absorbers are LLSs (17.3 < 
logA^(Hi) < 19.0), namely GRB 050908, GRB 060124, 
GRB 060607A, GRB 080310 - discussed in this paper ~, and 
GRB 090426. The GRB 080310 observations indicate that 
the line-of-sight does not cross any dense extended region 
with logiV(Hi) > 19. 

In our dataset, O i is in principle the most reliable metal- 
licity indicator, as it is a low-ionization species (in the same 
phase and location as the low-ionization gas) not strongly 
depleted onto dust. Moreover, ch arge exchange between O i 
and Hi (iField fc Steigma"nlll97lf ). caused by the similarity 
of their ionization potentials, assures a constant ratio of the 
two species, independently of ionization effects. However, 
this is only valid for a cloud in equilibrium, i.e., embedded 
in a co nstant radiat ion field, and perhaps not for log A^(H i) 
< 19 ()Viegaslll995D . Hydrogen and oxygen are ionized by 
the GRB afterglow at different rates in the case of neutral- 



hydrogen column-densities as low as logA^(Hi) ~ 18.7 (pa- 
per II). 

From the above, caution should be exerted when using 
Oi to estimate the metallicity of the GRB host galaxy in 
LLS or lower Hi column density systems, as this could un- 
derestimate the overall metallicity. Moreover, the peculiar 
abundances of the GRB 080310 absorber cannot be simpli- 
fied into one metallicity indicator (see Sect. 15. 4| ). 

Owing to the heavy saturation of the strong Lya transi- 
tion, even at this low column density, the transition cannot 
be split into separate components in the same way as the 
metal lines, preventing a component-by-component analy- 



5.2. Fell! origin and evidence for ongoing piioto-ionization 

The Feiii UV34 AAA1895, 1914,1926 triplet, arising from 
the Feili 17"^ excited levels (^83), has never been re- 
ported in a GRB afterglow until now. However, it has 
been previously observed in the spectra of extreme envi- 
ronments s uch as broad-absorption l ines BAL quasars and 
Carinae (iHartig fc Baldwin! Il986t IWampler et al.l 119951: 



'Baldwi n et al.l 



2000; H all et al 



19961 lLaoret al.1 Il997t [Johansson et all 



I2OO2D in emission or absorption. In partic- 



ular, these lines are detected in absorption in a rare class of 
iron l ow-ionization BAL quasars (FeLoBAL, iBecker et al.l 
119971: iHah fc Hutsemekerij |2003[) . In these environments 
Fell and Feiii excited levels are most likely populated via 
collisions, in a rather dense medium (particle density 10^" < 
n(Hi)/cm~^ < 10 ^^, column density logA^(Hi)-- 21.3, 
Ide Kool et alll2002[ ). However, the low iV(Hi) that we mea- 
sured for the GRB 080310 gas disfavours a collisional sce- 
naricP^. suggesting that a radiative process must be respon- 
sible for the population of the ^83 level. Moreover, the si- 
multaneous variability of Fe ill UV 34 and Fe 11 transitions 
favours a radiative process populating the Fe ill excited lev- 
els. 

In rj Carinae, the Feiii UV 34 triplet was observed in 
emission produced by the ion decaying from the term, 
being populated from the ground state via Lya pump- 
ing ([Johansson et al.ll2000l ). From this level, the ion has a 
high probability of decaying to the ''S3 level, and so in the 
presence of an increased Lya flux, this level can be popu- 
lated efficiently. However, we note that this Lya-pumping 
channel cannot play a relevant role in a GRB environ- 
ment, because both the GRB and pre-burst star-forming 
region cannot produce a significant Lya enhancement (pre- 
burst Lya photons cannot efficiently penetrate the absorb- 
ing cloud, while afterglow photons have little time to affect 
the line-of-sight Lya). 

In paper II, we present a novel explanation of the Fe ill 
UV 34 triplet along the GRB 080310 line-of-sight: in this en- 
vironment, the excited Feiii level ^83 is populated directly 
through the ionization of Fe 11. 

We note that the variability of Fe ill is stronger for com- 
ponent "b", while all the other features show a stronger 
variability in the total line profile (see Table [7]). This could 
suggest a higher rate of ionization of Fell into Feiii for 
component "b" , while the variability in the other compo- 
nents could be dominated by photo-excitation. Together 



® See Felll energy level diagram in I Johansson et aP (|2000l l. 
Even if the entire column would be compressed into a single 
parsec, the particle density would still only be ~ 2 cm~'^. 
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Fig. 6. Left: CLOUDY photo-ionization model used to derive the pre-burst physical conditions and metallicity in the 
GRB 080310 host-galaxy ISM, assuming the extragalactic UV background radiation field at the GRB host-galaxy redshift. 
In the main panel, column density is plotted against ionization parameter U for a range of metal lines. The observations 
are shown as large data points, colour-coded by ion. The best-fit log U (and hence the best-fit density) and [Fe/H] are 
derived by matching the observed Fciii/Feii ratio, as shown in the bottom panel. [C/H], [O/H], [Si/H], and [Cr/H] are 
determined by matching C ii, O i. Sill, and Cm at the best-fit log U. Right: Same as in the left panel but assuming the 
MW radiation field. 



with the more extreme values of [Si/Fe] and [C/Fe] in com- 
ponent "b", this could be another hint that this compo- 
nent is the closest to the GRB. A component-by-component 
GRB-cloud distance estimate with a photo-ionization and 
-excitation model could either confirm or refute this hy- 
pothesis, but the difficulty in distributing Hi between the 
individual components would render the results uncertain. 

In general, the observations of the GRB 080310 after- 
glow strongly indicate that we are witnessing the afterglow 
radiation gradually ionizing the surrounding medium. The 
pieces of evidence for such a photo-ionization in action are 
i) the strong variability of the Fell populations from all 
energy levels, in particular the ground-statcQ; the de- 
tection and variability of Fe ill transitions arising from the 
""Sa excited level, which are produced by ionization of Fc ii; 
Hi) the variability of the Cr ii ground-state population. The 
occurrence of photo-ionization in GRB 080310 is favoured 
by the low Hi content of the absorber, which does not effi- 
ciently shield the gas from the ionizing photons. However, 
it is mainly the high [Fc/H] that allows the Feiii UV34 
triplet to be detected (sec paper II). 



Although Sin and Fen have very simil ar ionization poten- 
tials, their photo-ionization cross-sections (|Vern6r et al.l [1993 . 
119961 ) are such that Sill is less ionized, which is consistent with 
the GRB 080310 observations. 



5.3. Nv non-variability 

The line profiles of the Nv doublet and the Fell fine- 
structure lines differ remarkably. In particular, the N v lines 
display two components around km s~^ and —50 km s~^, 
and clearly lacks the component at —20 km s~^. We do not 
detect any time variability in the Nv profiles (see Table [B]), 
which would have been expected had the lines bee n pro- 
duced within about 10 pc of the burst (IProchask a et al.l 
120081) . Thus, we exclude a circumburst origin for the ma- 
jority of Nv a bsorption i n GR B 080310. For QSO DLAs 
and sub-DLAs, iFox et al.l ()2009D showed that N v and C iv 
column densities are correlated and that their averaged b- 
values are similar, suggesting that the two ions arise in the 
same phase and have similar large-scale kinematics. Our 
analysis of the Nv line profiles and their non- variability 
favour the foreground host-galaxy ISM scenario for the pro- 
duction of Nv in the GRB 080310 afterglow. 

5.4. Relative abundances: 

ionization, supernova yields or dust destruction? 

Estimating relative abundances from ionic column densi- 
ties requires knowledge of ionization effects. Here, we first 
discuss relative ionic abundances (e.g., [Siii/Feii]; similar 
to the definition of the relative abundances of chemical ele- 
ments given in Sect.[T]) for ions that arc normally dominant 
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Table 8. Abundances towards GRB 080310 compared to other absorbers 



Ratio 


coinp. b 


comps. c+d 


total 


GRB DLAs 


QSO (sub-)DLAs 


QSO LLSs 


References 


[C/H] 


- 


- 


-1.3 + 0.2 






[-1.9, +0.2] 


[1] 


[0/H] 






< -0.8 










[Si/H] 


- 


- 


-1.2 + 0.2 


[-2.6, +0.7] 


[-2.6, +0.0] 


[-0.5, +0.0] 


[1],[2],[3],[4] 


[Cr/H] 


- 


- 


+0.7 + 0.2 




[-2.3, -0.8] 




[3] 


[Fe/H] 


- 


- 


+0.2 + 0.2 


[-3.0, +0.2] 


[-3.0, +0.0] 


[-0.6, -0.4] 


[2], [3], [5], [6], [7] 


[Ni/H] 


- 


- 


" < +0.95 




[-2.3, -0.6] 




[3] 


[Zn/H] 






" < +1.08 


[-1.8, +0.2] 


[-2.0, +0.0] 




[2], [4] 


[C/Fe] 


< -1.74 


< -0.91 


-1.5 + 0.2 




[-0.4, +2.0] 




[7] 


[O/Fe] 


< -1.09 


< -1.18 


< -0.8 




[-0.4, +1.2] 




[7] 


[Si/Fe] 


< -1.47 


< -0.72 


-1.4 + 0.2 


[-0.3, +1.2] 


[-0.1, +0.7] 




[2], [3], [7], [8] 


[Cr/Fe] 


" +0.56 ±0.19 


" < +0.79 


+0.5 + 0.2 




[-0.1, +0.5] 




[3], [8] 


[Ni/Fe] 


" < +0.45 


" < +0.70 


" < +0.52 




[-0.2, +0.3] 




[3] 


[Zn/Fe] 


" < +0.56 


" < +0.81 


" < +0.65 


[+0.0, +1.7] 


[+0.0, +1.0] 




[2], [7] 



Notes. Abundances with respect to either H or Fe estimated for the GRB 080310 absorber (see text for a description of the 
pre-burst ionization corrections). The ranges of abundances observed in GRB or QSO absorbers are listed on the right side of the 
table. All these abundances are visualized in Fig. [T] " Ionic relative abundances, i.e., not including the ionization effects of either 
the GRB or the ambient radiation field. ^ Does not include the possible - albeit limited - underestimate of Si and C in component 
"c" at epochs II and IV. 

References. [1] Q SO LLSs, iProchaska fc BurlesI (119991'): [21 GRB DLA s, ISchadv et all (l201ll): [31 fProchaska et~all (|200lh : [4] 
Wolfe et all moEi): [5] QSO sub-DLAs. iDessauges-Zavadskv eTal] ((20091 ): [6] QSO LLSs. iProchaskal (|1999l ): [7] QSO sub-DLAs, 



Peroux et al. 



20031 ): [8] rLedoux et~aD (|2002l 'l. 



in Hi clouds. Wc then refer to total element abundances 
after discussing the effects of ionization and deriving ion- 
ization corrections from CLOUDY modelling. 



The relative ionic abundances in the absorber have par- 
ticularly low [Siii/Feii] = -1.47 ±0.14 and [Cii/Feii] = 
-1.74 + 0.17 ratios in component "b" ([Siii/Feii] > -1.02 
for the total line profile and [Siii/Feii] > —0.72 in compo- 
nents "c+d") at epoch IL Despite the conservative lower 
limit to [Siii/Feii], we do not expect the ratio to be signif- 
icantly underestimated (less than 0.1 dex), given the con- 
straints on Sin at epochs I and III. so we can consider 
[Siii/Feii] ^ —1.0, ~ —0.7 for the total line profile and 
components "c+d", respectively. The observed [Siii/Feii] 
is difficult to explain since Sill and Fell have similar ion- 
ization potentials, as well as similar cosmic abundances. 
The typical value of [Siii/Hi] and high [Feii/Hi] (with 
a very low [Cii/Feii]) indicate that the low [Siii/Feii] is 
produced by an iron overabundance in the absorber rather 
than a silicon deficiency, which appears to be particularly 
marked in the peculiar component "b" . We can exclude any 
significant fraction of Si ii being hidden in higher-energy 
non-observable excited states, given the large energy gap 
bet ween the first odd terms a nd the high decay probabil- 
ity (jMartin fc Kaufmanlll994 ). Similarly, we can exclude a 
significant contribution of higher-excited states of oxygen 
because even if the non-detection limit for Oi** A1306 is 
the same as for O i* A1304 (same oscillator strength), O i** 
A1306 is typically significantly weaker than O i* A1304. 



5.4.1. GRB ionization effects 

The larger ov erall photo-ionization cros s-section of Fe ii rel- 
ative to Sin (jVerner et al.lll993l [19961 ) indicates that Fen 
is ionized more effectively than Si n by the GRB afterglow. 
Therefore, the overabundance of Fe II must have been even 
more pronounced before the onset of the burst. We found 
a pre-burst [Sin/Fen] = —1.1 for the total line profile, 
using a photo-ionization and -excitation modelling of the 
GRB radiation on the surrounding medium (paper II; see 
Sect. ID). 

5.4.2. Pre-burst ionization effects 

The absorber being a LLS, ionization effects due to the 
host-galaxy and extragalactic radiation fields may influence 
the ionic abundances. At the low-ionization end, neutral 
species such as Fe i and Si i are unexpected since they are 
easily ionized. In LLSs, Si i and Fe i are expected to be more 
than 10 00 times less abundan t than their singly ionized 
species ()Tumlinson et al.]|2011[ ). This is confirmed by our 
CLOUDY photo-ionization modeUing (logiV(Fei) < 11 for 
the best-fit ionization parameter). 

Regarding the intermediate-ionization species, while we 
can measure Fein to be at most as abundant as Fen, Sim 
cannot be constrained from the spectra because Si in A1206 
is blended with a strong C iv doublet associated with an in- 
tervening system. However, that Fe ii is more easily ionized 
than Sin reduces the possible amount of Sim that could 
contribute to the overall Si content. This is also confirmed 
by the higher Feiii/Feii than Sim/Sin (and Ciii/Cii) ex- 
pected from_aj]_eguilA^ modelling of 
LLSs (|Tumlinson et all 1201 11 ) . given the GRB 080310 Hi 
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column density. Thus, [(Siii+Siiii)/(Feii-|-Feiii)] should 
be low as well. Additional strong support for this is pro- 
vided by the low C ii content in component "b" . Indeed, 
C II can trace both Si ii and Si ill because C ii has a much 
higher ionization potential than Sin. Therefore, we can 
exclude ionization effects such as the presence of Sim as 
being the cause of the peculiar values of [Siii/Feii] and 
[Cii/Feii]. The results of our CLOUDY modeUing confirm 
the low abundance of Sim compared to Feiii (see Fig. [Bl). 

On the high-ionization side, Ovi, Nv, CiV, and Siiv 
have different component str uctures and app ear to arise in 
a separate ionization phase (|Fox et al.ll2008|) . In addition, 
including higher ionization species in the abundance calcu- 
lations would also require including Hil, which cannot be 
measured from the spectrum. Our CLOUDY modelling in- 
dicates that the pre-burst contribution of Fe iv, Si iv - and 
thus also the high-ionization species with higher ionization 
potentials such as C IV and O vi - to the total column den- 
sities is negligible. 

Finally, a photo-ionization model for a LLS with N(Hi) 
and redshift similar to GRB 080310 (i.e., photo-ionized by 
the same UV background) confirms that for solar relative 
abundances, [Sii i/Feii] is expect ed to be > for all ioniza- 
tion parameters ()Prochaskall999l ). Thus, the low [Siii/Feii] 
observed in GRB 080310 is due to peculiar abundances 
rather than ionization. These ionization arguments limit 
the total abundances of Si and C relative to Fe, to [Si/Fe] 
< [Siii/Feii] and [C/Fe] < [Cii/Feii]. We rely on these 
estimates for the relative abundances of individual compo- 
nents. 

For the total line profile, we can estimate the rela- 
tive abundances with our CLOUDY photo-ionization model 
(Sect.H]), which is based on the pre-burst column densities 
derived in paper II. To assess how reliable these results are, 
we first note that these column densities are somewhat close 
to what we found at epochs I and II, when the burst had less 
time to ionize the surrounding medium, and they do not 
depend dramatically on the GRB photo-ionization model 
assumption (see paper II). Second, the relative abundances 
derived from the CLOUDY model for the adopted extra- 
galactic background are similar to those found assuming 
the MW radiation field, strengthening our results. Thus, 
our estimates of the pre-burst ionization-corrected abun- 
dances are solid. 

The relative abundances estimated in the GRB 080310 
absorber - including the ionization correction discussed 
above - are listed in Table [SJ along with literature values 
for large samples of both QS0-DLAf[3 and sub-DLAs, 3 
LLSs, and 18 GRB absorbers. The same abundances are 
visualized in Fig. [71 suggesting that there is a heavy-metal 
enhancement in the absorbing gas, in particular Fe and Cr 
and possibly also Zn and Ni, given their high upper limits. 

5.4.3. Supernova yields 

We have shown that the GRB 080310 line-of-sight exhibits 
peculiar abundances, and that these cannot be explained 



We do not include the upper or lower limit s to the QSO-DLA 
abundances listed by iProchaska et al.l (|200lh . In particular, we 
exclude the outlier in the [Si/Fe] distribution, [Si/Fe] > —0.617 
(J0255-I-00 absorber at z = 3.9), because it has been estimated 
from a single Si line in the Lya forest (Sill A1193) in a spectral 
region where the continuum level is poorly constrained. 
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Fig. 7. Abundances relative to H (top panel) and Fe (bot- 
tom panel) towards GRB 080310 (diamonds and arrows), 
compared to other absorbers (rectangular regions), taken 
from Table[5](see text for a description of the pre-burst ion- 
ization corrections). The LLS [X/H] abundances are based 
on the observations of a few LLSs, and the corresponding 
[X/Fe] derived by combining the [X/H] ranges (the result- 
ing [X/Fe] ranges are probably overestimated). Thin arrows 
are ionic abundances for GRB 080310, i.e., uncorrected for 
ionization, observed at epoch II. The [C/Fe], [0/Fe], and 
[Si/Fe] estimated in GRB 080310 are much lower than pre- 
viously observed in DLAs, sub-DLAs, and LLSs. 



by ionization effects either owing to the afterglow or the 
extragalactic background radiation field. One possible ex- 
planation of the peculiar relative abundances, in partic- 
ular the [Si/Fe] ratio, is that nearby supernovae (SNe) 
may have enriched the surrounding medium with heavy 
elements. Peculiar chemical abundances have occasionally 
been observed in QSO-DLAs and explained with supernova 
(SN) yields (e.g., lAeafonova et al.l 2011 : Kobavashi et al.l 
1201 1[ ). In particular, Agafonpva et al.l ()2011 ) suggested that 
the low [Si/Fe] - -0.6^3 estimated in their z = 0.45 
absorber may be the remnant of the SN la explosion of 
a high-metallicity white dwarf (u sing th e SN la abun- 
dance patterns of |SMile_eialJ [2005: Ma zzah et al.l 120081: 
iTanaka et ani201lt) . 

The metal yields from core-collapse SNe and hyper- 
novae, including GRB- SNe, provide < [Si /Fe] < +2, de- 
pending on the model (jNomoto et al.l [20091) . failing to ex- 
plain any sub-solar [Si/Fe] ratio. Mod els of Type la SNe 
predict [Si/Fe] ratios dow n to ~ -0.2 ([Maeda et al.ll2010t) 
or as low as to ~ —0.5 (jTanaka et al.l 12011 ) 



for Type la 

SN chemical enrichment, together with slightly super-solar 
[Cr/Fe] and [Ni/Fe] and low [Zn/Fe] ratios. 



This Si estimate should be taken with caution, because it 
is derived from Sir assuming that Sil/Sill = Mgl/Mgll, which 
could lead to a significant Si underestim ate, see photo-ionization 
models in, e.g.. iTumlinson et al.l (|201lh . 
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These predictions are partially consistent with the rel- 
ative abundances observed in components "c-|-d" of the 
GRB 080310 absorber, suggesting that SNe la may have 
enriched the ISM with heavy metals along the line-of-sight. 

However, the extremely low [Si/Fe], [0/Fe], and [C/Fc] 
measured in component "b", as well as the pre-burst val- 
ues inferred for the total line profile, cannot be repro- 
duce d only with SN la yie lds and require further explana- 
tion. lOhkubo et all (|2006D found that a low [Si/Fe] 2.5 

could be expected for the chemical enrichment by the core- 
collapse of very massive Pop III stars (500-1000 Mq). 
However, these models also predict very high Ni and Zn 
contents ([Ni/Fe] ~ -1-1.0, [Zn/Fe] ^ +1.6) and a strong 
depletion of chromium ([Cr/Fe] ~ —1-0), all features that 
we do not observe. Exotic processes such as ONcMg core 
explosions or accretion-induced white-dwarf collapses could 
qualitatively explain these peculiar abundances, but they 
could only contribute marginally to the ISM abundances. 

In general, the low [Si/Fe] observed in the GRB 080310 
absorber also requires no a-element contributions from 
core-collapse SNe, indicating that there has been a small 
amount of recent star formation in the host galaxy. Thus, 
a long G RB (likely aris i ng from a massive sta r progeni- 
tor, e.g., IWooslevlTTool iHiorth fc Bloomll201lh indicates 
episodic massive-star formation in the host galaxy. 

The lack of recent star formation along the line-of-sight 
is perhaps another indication - after the low A^(H i) - that 
GRB 080310 is not located within the densest parts of its 
host galaxy, given that GRB hosts are typically actively 
star-forming (e.g., ISavaglio et al.l [20091 ). The evidence of 
some GRB absorbers - those with almost featureless con- 
tinua - that are not deeply embedded in the densest regions 
of their host galaxies was discussed in lDe Cia et al.l (|201lt) . 
However, several, strong absorption lines are detected in the 
case of GRB 080310. This could suggest that the LLS en- 
vironment of GRB 080310 differs intrinsically from those of 
most GRB-DLAs, corresponding to possibly a neutral-gas 
poor and not highly star-forming host galaxy, and that the 
measurements are not simply due to a particular location 
of the burst within its host galaxy. 

5.4.4. Dust destruction 

One intriguing possibility for explaining the Fe and Cr 
enhancements in component "b" could be dust destruc- 
tion, possibly induced by the GRB itself. From the Voigt- 
profile modelling. Fell and Cm show stronger absorption 
in component "b", suggesting that these ions are mostly 
co-spatial. Moreover, Fe and Cr are the most depleted el- 
ements in the dust grains of QSO- and GRB-DLAs (e.g., 
ISavaglio et al.l[2C)03t) . The destruction of silicate grains can 
release Fe and recycle it into the ISM out to at least 3 pc 
from the GRB, increasing iV(Fe) by one order of magnitude 
and not a ffecting the abundance of less depleted elements 
such as O (jPerna &: Lazzatill2002[ ). In addition, Fe and FeaC 
are the first and sometimes the onl y grains to conden sate 
at high temperature (T - 1000 K. iLcwis fc Nevll 19791) . 

If dust destruction recycles the metals into the ISM, 
the relative abundances of the gas no longer suffer from 
dust depletion and represent the intrinsic abundances. This 
alone does not explain the peculiar Fe and Cr overabun- 
dances that we observe. However, this only applies when 
the gas and dust evolve together, which is not necessar- 
ily the case. In particular, radiation pressure can cause the 



dust grains to drift away from the inner parts of H ii regions 
(|Draindl20lTh . Large dust grains can also be driven out to 
at least beyond 0.5 pc into the ISM by stellar- wind shocks 
(Ivan Marie et al.ll201ll) . Both of these processes can there- 
fore change the local dust/gas ratio within star- forming re- 
gions. The iron overabundance in component "b" may be 
the result of the destruction of Fe-rich dust grains, such as 
iron silicates or solid Fe, while the low C abundance sug- 
gests that this destroyed dust cloud has a low carbonaceous 
content. 

The GRB radiation can quickly destroy dust grains out 
to a distance of ^100 pc from the burst (jFruchter et al.l 
I2OOII ). Alternatively, a GRB-unrelated process such as 
shocks (e.g., in S N ejecta) could be responsible for the dus t 
destruction (e.g., iJones fc NuthI 120111: iJones et al.l I1996D . 
This has never been reported in QSO absorbers, but QSO 
lines-of-sight typically probe the outer regions of the in- 
tervening galaxies, where a low dust content and SN rate 
are expected. In conclusion, the destruction of drifted Fe- 
rich dust grains, possibly due to the burst or a nearby SN, 
could explain the iron overabundance in component "b" , 
but further modelling is clearly needed to investigate the 
feasibility of this scenario. 

In any case, a large amount of dust in the absorber can 
be excluded, because if a significant amount of iron would 
be depleted into dust then the intrinsic [Si/Fc] would have 
to be even lower. 



6. Conclusions 

Our time-series of high- resolution VLT/UVES spectra of 
the GRB 080310 afterglow has revealed the unique features 
of the ISM of the GRB host galaxy. We have reported the 
detection of several resonance absorption lines commonly 
observed in the ISM, as well as Fell fine-structure lines 
sometimes associated with GRB afterglows. In the case 
of GRB 080310, both the Fell ground-state and the fine- 
structure lines vary with time. We decomposed the complex 
spectral-line profiles of several ions and modelled them al- 
together using a four-component Voigt profile, resulting in 
column-density determinations for each component. These 
components might be associated with different clouds along 
the line-of-sight within the GRB host galaxy. Interestingly, 
a low logA^(Hi) = 18.7 ± 0.1 is derived from the Lya ab- 
sorption. 

We also detected the Feiii UV34 AAA1895, 1914,1926 
line triplet in absorption. These transitions arise from the 
^Ss energy level of Fe ill, which is observed for the first time 
in a GRB afterglow spectrum. The Fe il and Fe ill time vari- 
ability, of both the ground-state and excited levels, is clear 
evidence that we are witnessing the ISM being gradually 
photo-ionized by the GRB afterglow radiation. 

Through the analysis of chemical abundances mea- 
sured in the GRB 080310 absorber and CLOUDY photo- 
ionization modelling, we inferred pre-burst ionization- 
corrected ratios of [C/Fc] = -1.5 ± 0.2, [0/Fe] < -1.0, 
and [Si/Fe] = -1.4 ± 0.2 (where [Fe/H] = +0.2 ± 0.2 and 
[Cr/H] = -f0.7±0.2), while, typically, < [Si/Fe] < +0.7 is 
observed in QSO and/or GRB absorbers. Furthermore, we 
observed an even more extreme iron overabundance ([Si/Fe] 
< —1.47 and [C/Fe] < —1.74) in the gas associated with 
component "b" of the absorption profile. Such a low [Si/Fe] 
ratio has never been observed before in cither QSO- or 
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GRB-DLAs and cannot easily be explained by current mod- 
els of SN la and Pop III SN chemical yields. A potential 
explanation might be provided by the destruction of iron- 
rich dust grains, thereby recycling heavy elements into the 
gas phase. The dust could be destroyed by the GRB itself 
or GRB-unrelated processes such as SNe shock waves. 

The high iron column-density measured in the gas phase 
generally suggests a low dust content in the absorber. The 
strong overabundance of iron compared to silicon and car- 
bon also suggests that there has been negligible recent star 
formation along the line-of-sight. The occurrence of the 
GRB then indicates that there has been episodic massive- 
star formation in the GRB region. 

Acknowledgements, We thank Peter Laurson and Keiichi Macda for 
insightful discussions and Jason Prochaska and Sandra Savaglio for 
useful comments. We also thank an anonymous referee for a care- 
ful and constructive report. ADC acknowledges support from the 
ESO Director General Discretionary Fund 2009 and 2010, and the 
University of Iceland Research Fund. PJ acknowledges support by a 
Marie Curie European Re-integration Grant within the 7th European 
Community Framework Program and a Grant of Excellence from the 
Icelandic Research Fund. The Dark Cosmology Centre is funded by 
the Danish National Research Foundation. We acknowledge the care- 
ful assistance of the VLT observers, in particular Claudio Melo and 
Dominique Naef. 



References 

Agafonova, I. I., Molaro, P., Levshakov, S. A., & Hou, J. L. 2011, 
A&A, 529, A28 

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 
47, 481 

Baldwin, J. A., Ferland, G. J., Korista, K. T., et al. 1996, ApJ, 461, 
664 

Ballester, P., ModigUani, A., Boitquin, O., et al. 2000, The Messenger, 
101, 31 

Becker, R. H., Gregg, M. D., Hook, I. M., et al. 1997, ApJ, 479, L93 
Christensen, L., Fynbo, J. P. U., Prochaska, J. X., et al. 2011, ApJ, 
727, 73 

Cummings, J. R., Baumgartner, W. H., Beardmore, A. P., et al. 2008, 

GRB Coordinates Network, 7382 
De Cia, A., Starling, R. L. C, Wiersema, K., et al. 2011, MNRAS, 

418, 129 

de Kool, M., Becker, R. H., Arav, N., Gregg, M. D., & White, R. L. 

2002, ApJ, 570, 514 
Dekker, H., D'Odorico, S., Kaufer, A., Delabre, B., & Kotzlowski, 

H. 2000, in Society of Photo-Optical Instrumentation Engineers 

(SPIE) Conference Series, Vol. 4008, , 534-545 
Dessauges-Zavadsky, M., Ellison, S. L., & Murphy, M. T. 2009, 

MNRAS, 396, L61 
Draine, B. T. 2011, ApJ, 732, 100 

Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 
761 

Field, G. B. & Steigman, G. 1971, ApJ, 166, 59 

Fiore, F., D'Elia, V., Lazzati, D., et al. 2005, ApJ, 624, 853 

Fontana, A. & Ballester, P. 1995, The Messenger, 80, 37 

Fox, A. J., Ledoux, C, Vreeswijk, P. M., Smctte, A., & Jaunsen, A. O. 

2008, A&A, 491, 189 
Fox, A. J., Prochaska, J. X., Ledoux, C, et al. 2009, A&A, 503, 731 
Fox, A. J., Wakker, B. P., Savage, B. D., et al. 2005, ApJ, 630, 332 
Fruchter, A., Krolik, J. H., & Rhoads, J. E. 2001, ApJ, 563, 597 
Fynbo, J. P. U., Jakobsson, P., Prochaska, J. X., et al. 2009, ApJS, 

185, 526 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Hall, P. B., Anderson, S. F., Strauss, M. A., et al. 2002, ApJS, 141, 
267 

Hall, P. B. & Hutsemckers, D. 2003, in Astronomical Society of 
the Pacific Conference Scries, Vol. 290, Active Galactic Nuclei: 
From Central Engine to Host Galaxy, ed. S. Collin, F. Combes, 
& I. Shlosman, 209 

Hartig, G. F. & Baldwin, J. A. 1986, ApJ, 302, 64 

Hjorth, J. & Bloom, J. S. 2011, to appear in "Gamma-Ray Bursts", 
eds. C. Kouveliotou, R. A. M. J. Wijers, S. E. Woosley, Cambridge 
University Press (ArXiv: 1104.2274) 



Jakobsson, P., Fynbo, J. P. U., Ledoux, C, et al. 2006, A&A, 460, 
L13 

Johansson, S., Zethson, T., Hartman, H., ct al. 2000, A&A, 361, 977 
Jones, A. P. & Nuth, J. A. 2011, A&A, 530, A44 

Jones, A. P., Tielens, A. G. G. M., & Hollenbach, D. J. 1996, ApJ, 
469, 740 

Kobayashi, C, Tominaga, N., & Nomoto, K. 2011, ApJ, 730, L14 
Laor, A., Jannuzi, B. T., Green, R. F., & Boroson, T. A. 1997, ApJ, 
489, 656 

Ledoux, C, Bergeron, J., & Petitjean, P. 2002, A&A, 385, 802 
Ledoux, C, Vreeswijk, P. M., Smette, A., et al. 2009, A&A, 506, 661 
Lewis, J. S. & Ney, E. P. 1979, ApJ, 234, 154 

Lodders, K., Palme, H., & Gail, H.-P. 2009, in "Landolt-Bornstein 
- Group VI Astronomy and Astrophysics Numerical Data and 
Functional Relationships in Science and Technology Volume, ed. 
J. E. Triimper, 44 
Maeda, K., Ropke, F. K., Fink, M., et al. 2010, ApJ, 712, 624 
Martin, W. C. & Kaufman, V. 1994, VizieR Online Data Catalog, 
6074, 

Mazzali, P. A., Sauer, D. N., Pastorello, A., Bcnetti, S., & Hillcbrandt, 

W. 2008, MNRAS, 386, 1897 
Meszaros, P. 2006, Reports on Progress in Physics, 69, 2259 
Nomoto, K., Wanajo, S., Kamiya, Y., Tominaga, N., & Umeda, H. 

2009, in lAU Symposium, Vol. 254, , 355-368 
Ohkubo, T., Umeda, H., Maeda, K., et al. 2006, ApJ, 645, 1352 
Pcrna, R. & Lazzati, D. 2002, ApJ, 580, 261 
Pcrna, R. & Loeb, A. 1998, ApJ, 501, 467 

Peroux, C, Dessauges-Zavadsky, M., D'Odorico, S., Kim, T.-S., & 

McMahon, R. G. 2003, MNRAS, 345, 480 
Piran, T. 2004, Reviews of Modern Physics, 76, 1143 
Prochaska, J. X. 1999, ApJ, 511, L71 
Prochaska, J. X. & Buries, S. M. 1999, AJ, 117, 1957 
Prochaska, J. X., Chen, H., Dessauges-Zavadsky, M., & Bloom, J. S. 

2007, ApJ, 666, 267 
Prochaska, J. X., Dessauges-Zavadsky, M., Ramirez- Ruiz, E., & Chen, 

H.-W. 2008, ApJ, 685, 344 
Prochaska, J. X., Foley, R., Tran, H., Bloom, J. S., & Chen, H. 2006, 

GRB Coordinates Network, 4593 
Prochaska, J. X., Wolfe, A. M., Tytler, D., et al. 2001, ApJS, 137, 21 
Ranch, M. 1998, ARA&A, 36, 267 

Rothman, L. S., Gordon, I. E., Barbe, A., et al. 2009, 

J. Quant. Spec. Radiat. Transf., 110, 533 
Sargent, W. L. W., Steidel, C. C, & Boksenberg, A. 1989, ApJS, 69, 

703 

Savaglio, S., Fall, S. M., & Fiore, F. 2003, ApJ, 585, 638 
Savaglio, S., Glazebrook, K., & Le Borgne, D. 2009, ApJ, 691, 182 
Schady, P., Savaglio, S., Kriihler, T., Greiner, J., & Rau, A. 2011, 
A&A, 525, A113 

Smette, A., Sana, H., & Horst, H. 2010, Highlights of Astronomy, 15, 
533 

Stehle, M., Mazzah, P. A., Bcnetti, S., & Hillebrandt, W. 2005, 

MNRAS, 360, 1231 
Tanaka, M., MazzaU, P. A., Stanishev, V., et al. 2011, MNRAS, 410, 

1725 

Thone, C. C, Campana, S., Lazzati, D., et al. 2011, ArXiv c-prints 
Tumlinson, J., Werk, J. K., Thorn, C, et al. 2011, ApJ, 733, 111 
van Marie, A. J., Meliani, Z., Keppens, R., & Decin, L. 2011, ApJ, 
734, L26 

Vergani, S. D., Petitjean, P., Ledoux, C, ct al. 2009, A&A, 503, 771 
Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, 
ApJ, 465, 487 

Verner, D. A., Yakovlev, D. G., Band, I. M., & Trzhaskovskaya, M. B. 

1993, Atomic Data and Nuclear Data Tables, 55, 233 
Viegas, S. M. 1995, MNRAS, 276, 268 

Vogt, S. S., Allen, S. L., Bigelow, B. C, et al. 1994, in Presented 
at the Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference, Vol. 2198, SPIE Conference Series, ed. D. L. Crawford 
& E. R. Craine, 362 

Vreeswijk, P. M., Ellison, S. L., Ledoux, C, et al. 2004, A&A, 419, 
927 

Vreeswijk, P. M., Jakobsson, P., Jaunsen, A. O., Oslo, U., & Ledoux, 

C. 2008, GRB Coordinates Network, 7391 
Vreeswijk, P. M., Ledoux, C, Raassen, A. J. J., et al. 2012, A&A, 

000, accepted 

Vreeswijk, P. M., Ledoux, C, Smette, A., et al. 2007, A&A, 468, 83 
Wampler, E. J., Chugai, N. N., & Petitjean, P. 1995, ApJ, 443, 586 
Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861 
Wolfe, A. M., Turnshek, D. A., Smith, H. E., & Cohen, R. D. 1986, 



16 



A. De Cia et al.: Rapid-response mode VLT/UVES spectroscopy of super iron-rich gas exposed to GRB 080310 



ApJS, 61, 249 
Woosley, S. E. 1993, ApJ, 405, 273 



Appendix A: 



A. De Cia et al.: Rapid-response mode VLT/UVES spectroscopy of super iron-rich gas exposed to GRB 080310 
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Fig. A.l. Voigt-profilc fit of absorption lines observed at epoch II. Tire fit decomposition was derived from the highest 
S/N spectrum (obtained at epoch IV). 
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Fig. A. 2. Voigt-profile fit of absorption lines observed at epoch I. The fit decomposition was derived from the best S/N 
spectrum (obtained at epoch IV). For the most poorly resolved lines, e.g., Feiii A1122, we conservatively adopt upper 
limits, see Table E) 



Table A.l. Intervening systems: identified absorption lines 



z 


I'host 


Transitions 




VK^^est (A) 










A (A) 


Lya 


C IV 1548 


Mgii 2796 


Sin 1526 


Fe II 2382 


2.4200 


650 


Lya, Civ 1548, 1550^ 
Fe II 2382 


1.04 ±0.01 


0.05 ±0.01 


< 0.05 


< 0.02 


0.05 ±0.01 


2.4113 


1411 


Lya, Civ 1548, 1550^= 


0.70 ±0.01 


0.13 ±0.02 


c 


< 0.03 


< 0.02 


2.2786 


13018 


Lya, Cll 1334, Cll* 1335, 
Civ 1548,1550, 
Mgii 2796, 2803, 
Sin 1260^ Sim 1206, 
Si IV 1393, 1402, 
AIii 1670, Aim 1854, 1862 


1.62 ±0.01 (« 18. 7^*) 


1.03 ±0.02 


0.07 ±0.02 


< 0.02 


< 0.02 


2.1702 


22499 


Lya, CiV 1548, 1550 


0.61 ±0.01 


0.07 ±0.01 


c 


< 0.02 


< 0.02 


2.0685 


31395 


Lya, Civ 1548, 1550 


1.22 ±0.02 18.4^*) 


0.21 ±0.01 




< 0.02 


< 0.02 


1.6711 


66155 


Lya, Cll 1334, Cll* 1335, 
Civ 1548^1550, 
Mgii 2796, 2803, 
Alii 1670, Aim 1854, 1862, 
Sin 1526, Si IV 1393, 1402, 
Fell 1608, 2344, 2382, 2600 


< 0.96 


c 


0.42 ±0.01 


0.13 ±0.01 


0.04 ±0.01 


1.1788 


109216 


Sin 1526, 
Mgii 2796, 2803 




< 0.18 


0.05 ±0.01 


0.24 ±0.04 


< 0.04 



Notes. " Blue-shifted velocity with respect to the z = 2.42743 host-galaxy system. ^ Rest W and la error estimate, or 3(7 upper 
limit for non-detections, measured in the highest S/N spectrum for each line. The assumed 2xFWHM apertures correspond to 
a b-value of 20 km s~^, which enabled us to cover each line (except for Civ 1548 at z = 2.41, 2.27, 2.06, 1.67, where an aperture 
corresponding to fe = 100 km s~^ was necessary and fe = 50 km s~^ for Mgll 2796 at 2 = 1.67 and C iv 1548 at z = 1.17). Blended 
with another line. log7V(Hl) estimate from the W, assuming a damped regime {N-^ ^[cm~^] ~ 1.87 x 10^**M^^[A]) 
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Table A. 2. Intervening system at z = 1.6711: column densities of Alii and Aim 



Ion 


log N ± 


a 




a 




Alii 


11.33 ±0.26 


12.26 ± 0.07 


Aim 


12.28 ±0.19 


13.01 ± 0.04 


z 


1.671110 ±0.000007 


1.670923 ± 0.000029 


b 


8±5 km s"^ 


9 ± 1 km s"^ 



Notes. " Logarithm of the ion column density and corresponding Icr error estimate, for components a and /3. 




-50 50 

Velocity (km s'') 

Fig. A. 3. Voigt-profile fit of Alii and Aim - associated with an intervening system - and the host-galaxy Ol, observed 
at epoch IV, assuming the best-fitting log7V(Oi,"b")= 13.50 value. The Alii A1670 l ine at z = 1.6711 is blended with 
O I A 1302 at the GRB host-galaxy redshift. The O i line estimate is discussed in Sect. 13.2.31 
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